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ABSTRACT 


Beginning with linearized forms of the vorticity 
equation and the first law of thermodynamics applied to a 
Simple baroclinic model including friction, an analytic 
solution is obtained which describes the time variation 
of the relative phase angle between the temperature and 
pressure waves in the atmosphere. It is found that for 
unstable waves the atmosphere tends toward a state where 
the temperature wave lags the pressure wave. The time 
variation of the amplitude and relative phase depends on 
the initial value of the phase difference, but the ultimate 
angle by which the temperature field lags the pressure 
field depends only on the drag coefficient, certain atmos- 
pheric parameters, and the wavelength of the waves. Figures 
are included to show the variation with time and initial 
phase difference of the amplitude and the relative phase. 
Also shown are the effects on the amplitude and terminal 
value of the phase difference due to variations in the drag 
coefficient and other parameters. 

The writer wishes to express his appreciation to 
Professor George J. Haltiner of the U. S. Naval Postgraduate 
School for his guidance, contributions and encouragement 


in this work. 


alial 


Section 


TABLE OF CONTENTS 


Title 
Introduction 
Development 
Discussion of Results 


Bibliography 


Seals 


Page 


© 
22 








Figure 


ie 


LIST OF ILLUSTRATIONS 


Effects of variations 
the terminal value of 


Effects of variations 
the amplitude 


Time variation in the 


Effects of variations 
amplitude 


Effects of variations 


in the drag coefficient on 
the relative phase lag 


in the drag coefficient on 


amplitude of stable waves 


in the thermal wind on the 


in the static stability 


parameter on the amplitude 


Effects of variations 
amplitude 


in the latitude on the 


Time variation and dependence on « of the relative 


phase lag for a 5,000 


km wave 


Time variation and dependence on « of the amplitude 


for a 5,000 km wave 


iv 


Page 
14 


Le 


16 


My 


18 


1S) 


20 


Pall, 








Candace 


the 


time (dynes cm7 


the 
the 
the 
the 
the 
the 


the 


TABLE OF SYMBOLS 


total derivative 5 pressure with respect to 
sec™ 


atmospheric pressure (dynes cm7* 


) 
relative worticity. a> av —ad 
wind velocity vector: Y=yl+vf 
thermal wind vector 
coriolis parameter (sec71) 

finite difference of the parameter x 


average value of the coriolis parameter. and 


is equal to the value at latitude 45N (sec71l) 


the 
the 
the 
the 


the 


(cm) 


the 


the 


gradient operator: V=L ye +I3F¢7 
static stability parameter (cm? dynes~*) 
thickness of an atmospheric layer (cm) 
potential temperature (°K) | 


vertical height coordinate; Z=O at the surface 


magnitude of the vector {YW (cm sec7+) 


vertical component of the acceleration due to 


gravity (cm sec 


the 
the 
the 
the 
the 
the 
the 


surface drag coefficient (dyne sec em7?) 
stream function: YW= oh 

wave number: “=277/L 

wavelength (cm) 

phase speed of the waves 


Rossby parameter: fat 


x and y unit vector respectively 





~~ 
a 
> 
= 
— << = a 
= 
- /  —— 
_ @ 
- a ae 
7 = - 
a _ -_ =z 
-_> > «a 
ss — 





the -/-7’ 


the Rossby wave speed 


the amplitude of the pressure wave at t=0 


the amplitude of the temperature wave at t=0 


the relative phase between the temperature and 


pressure 
the time 
the time 
the time 


the time 
wave 


the basic zonal thermal wind (cm sec 


the basic zonal wind (em sec 


time 


wave at t=0 

varying amplitude of the pressure wave 
varying amplitude of the temperature wave 
varying phase angle of the pressure wave 
varying phase angle of the temperature 

7) 

-1) 


vi 






i. Introduction 

In 1960 Wiin-Nielsen [5] published the results of 
an investigation of the mechanisms in the atmosphere 
whereby the conversion of potential to kinetic energy 
is possible. He stated that for this conversion to occur 
it was necessary that the temperature wave lag behind the 
pressure wave. Using a simple baroclinic model with no 
friction he obtained an analytic solution for the relative 
phase between the temperature and pressure fields. Based 
on this solution he indicated how the relative phase angle 
changed with time for various initial values and how the 
early tendency of the pressure wave to amplify was affected 
by the initial value of the relative phase angle. 

The purpose of this work is to include friction in the 
simple baroclinic model and compare the results to the 
non-friction case. Also the values of the initial relative 
phase angle, wavelength, drag coefficient, Coriolis para- 
meter, thermal wind, and static stability were varied to 
investigate their effects on the amplification of the 
pressure wave, and the time variations including the limit- 
ing value of the relative phase between the temperature 


and pressure waves. 








en Development. 

The atmosphere is represented by a simple baroclinic 
model. The 750- and 250-mb levels are used as data levels, 
and the 500-mb level is taken as the level of interest. 
This gives a 4-layer model where @=0 at the upper boundary, 
P=0, (level zero); and @)=4), , the frictionally induced value 
at level four, the lower boundary of the atmosphere. It is 
also assumed that the atmospheric motions are adiabatic. A 
zonal basic wind field which is independent of latitude and 
a linear function of pressure is utilized. 

The equations used are the vorticity equation in the 


form 


23+: V(it-f, $2=0, (1) 


and an expanded form of the first law of thermodynamics for 


adiabatic motion given by Haltiner [2] 


a +Y-Vh =ANPaoa=o. 


Here o> is a measure of the static stability, 
oan L dz 26" 
ee age 
and ae has been replaced by —h/AP. It will be assumed 
that the basic flow pattern remains constant and that the 
perturbations grow, but no attempt will be made to account 


for the form of the necessary energy input for this to occur 


notwithstanding the dissipating effect of the friction. The 








vorticity equation will be applied at 25C mb, level one, 

and at 750 mb, level three. The thermal equation is applied 
at 500 mb, level two. The wind at level i is denoted Me, 
except that at level two the wind is taken to be V=4(Y.+Y) 
The thermal wind for the layer of thickness AP centered 
about level two is given by Y=F(Y-V)- 

For simplicity in the equations it was assumed that 
the frictional stress is directly proportional to the 
surface wind. The geostrophic wind approximation leads to 
the result a,=- 20, % Here C, is a drag coefficient, 
a5 is geostrophic relative vorticity, g is gravity, and f 
is the Coriolis parameter. 

Applying the vorticity equation at levels one and 


three yields: 


SL + Y WS +f) = fe $s; 





=f Gieee 
2 +4 Vie tf) = hn SS 
and expressing VY, and VW in terms of V and \% yields 


for level one the result 


(Th) +(V 4) WIT +6) = fe Sp (3) 


and for level three 


£-(f-5)+(V-W VT eh SF. ow 


Here, the thermal equation is applied at level two to obtain 


@ Seer: +VY-Wh) (Cea?) 











and the linear extrapolation YW +(BA)_AP is used to obtain 


a), 


BStV VT) +YyeVG =e (F-2). ©) 


Subtracting equation (4) from equation (3) produces 


525 tVe IO V-VS- eB + V-h)+ 
ere to) 


TGF OP 


To linearize these last two equations assume that YG and 


Adding equations (3) and (4) gives 


V are given by 


Pn) ee, and Y= +u)e +s, 

where the U and U. represent the basic zonal flow which is 

a function of the pressure only, and the small u's and v's 
represent the perturbation quantities which are independent 
of latitude, the y-axis. Also take h to be a basic part 
plus a perturbation, i.e., h=AW thy). Now by neglecting 
terms involving products of perturbation quantities and 
replacing v by vey, where (# is a stream function given 

by pat Z23 and V7 by Yate, where YW is a stream function 
given by weth, we obtain two third order linear partial 
differential equations in W and ce 


2 (Husk +l 2% +f Bee 3h +5 $e - 


— SUge NG 2 Se )=0 (7) 
35 (SX) +Usk +U, FE Ap (= “22% )=0 (8) 


where. (32 of. yO N= —fm'ICo/#FOP, and S=thn/2I6(Op)- 





Assume solutions of equations (7) and (8) of the form 


ty (x-Ct) Lu(x-Ct) 
W=Ke and Vr=Yoe 
= e4 
where Ua and VY, are complex amplitudes, 4/= UL is the wave 
number for wavelength equal to L , and C is the phase 
speed. Substituting these two solutions into the differ- 


ential equations yields two homegeneous linear equations 


in Y and Y, 


(C-Ur 4-22) % -(-¢B) Yo= 0 (9) 
(U(e-+iB)]¥ 4 +(C- ee +32 - ot -< é Aly <0 (enep 


In order to have non-trivial solutions to this system of 
equations it is required that the determinant of the 
coefficients of ¥ and Ye, ve equal to zero. Setting the 
determinant equal to zero and defining Cp=//- ay the 
Rossby wave speed, yields a quadratic equation in (C-Cr) 


Solving this equation for (€-Cp) gives 


C-Ge=F tiE ty [HuG' i 
F=BS/2auers) — £= (ND uss) 


G=[H4(2¢- a3) +28 (/-£) (aby 


H= vl -jx)-2 (5+ +a Az(i+®) 


By comparing equation (11) to the result given by 


Wiin-Nielsen [5] it can be seen that the effect of including 





friction in the model has been the addition of those terms 
involving A as a factor to the frequency equation. If 

the drag coefficient, Cp , is set equal to zero, A vanishes 
in equation (11) and the result is identical to Wiin- 
Nielsen's equation. 

Equation (11) gives C as a complex number where the real 
part determines the speed of propogation of the perturba- 
tions, and the imaginary part when multiplied by the wave 
number gives the amplification for the wave. Taking first 
the plus sign and then the minus sign with the radical gives 
two solutions which are designated as C,and C_ respectively. 


The complete solutions for ~and Ware then given by: 


VW a, CMS 2) 4+ DeH(x-c.t) (a2)) 


Ga GCMs rh CHO-c2) (23) 


Next take as anitial conditions at time t=o 


A 4 

a £4 X = Zz +r oX 
Wx CHR and Ue=Ap CLuUxr ed) 
where IM and An are real numbers and ec is a relative 
phase angle which is positive when the temperature field 
lags the pressure field. Then by considering equations (12) 
and (13) at time zero it follows immediately that 

ea) Lo 

Wty =A ad K.th=a4e 
By utilizing equation (9) for C+ and C_ we obtain two more 
equations for Yor, Yo-y Vros, ane V4,. giving a system of four 


equations and four unknowns. Solving the system gives 


Yh =(KAr Oo UA' Xe -X-) Vos 
YEKAC A (x-X.) Vo 
where K=lf> -¢ 2 2a 

X+=Ci-Cra-Z a and X= CSAGEE Se - 


Then by substituting the expressions for Yor, “-, Ve 


Y, (14) 
=y (25) 


4 xh 


and YW,_ into equations (12) and (13) the complete solutions 
for Wand YF are obtained. 

In order to study unstable baroclinic waves in some 
detail, it is desirable to have the expressions for ~and 
Vr in slightly different forms. Since C; and C_ are complex 
they can be expressed as C,;=YtiW and C_=M+iN. Then the 


expression for Wbecomes 


2 Awt_ L£H(x-ye) € LM x=. 
A ene Ye eee ae 
Now define Bt, B-, ps, and f- by 


VA apie 1m Wie Gye eS 


Then 


ss LF _uM(x-Ye) LR. 5 tMx—t) 
YU Be Cc +13-C ce 5 


which after some manipulation can be written as 
Y yRie £LAax—M(Cet PF) t +P] 


Here {/3is given by 


B*= Be +B. + 213A Com l(,-B)-4] 127) 





<pis given by 

Pee ig (34 sin(e-E)+ Boo (+E) (18) 
Been (Py pe Pre) 7 

Wis defined as Ha ZirlenG t where Y= Gea" and 


Q=arctan (G/H); F, G, H, and A are defined immediately 
following equation (11). The treatment for Wis similar 
to the above procedure for Y% Define Bre, Br, P, and 2 by 


Wwe LA Bs 
% en =/37,,€ °° and Ww, coe Br. © ae (19) 


7 


and then obtain 
Y%=Be LLAX—M Ont) L + Pr] 


with B,-given by 


(rie Ep +f, +2 [Sp [Sp- aM AH) iz) 


and given by 


— Br, se oe -$)+Bp Sl. + £) (25 


Br+ com R- ¥) + Br cor(4_+) 


We now have obtained expressions which give the variations 
of Wand yeas simple sinusoidal waves having amplitudes and 
phase angles which are functions of time. 

We now want to determine how the relative phase angle 
between the thermal and pressure waves, (A— > ) changes 
with time. Considering the exponential nature of Bs B-, Crys 
and [37., it can be shown from equations (18) and (21) that 


as time approaches infinity, op approaches (g-F) and Pf, 


approaches (#,-#). Therefore, in the limiting case 
(9--P) =(4,.-&). Now by suitable combining equations 
(14), (15), (16), and (19), it can be shown that 


Lan ( BA, a F) = 
§— GAC ste 2 1Y-N) 2 (BE Velcon 2 


~ UBM +P lees. 2)—2A (MU db ip se B) 7) 


Equation (22) was obtained by considering only the 
amplifying part of the wave. Unstable baroclinic waves 
tend toward a state such that the limiting phase angle 
between the thermal and pressure waves, (F.- FA), is given 
by He [ ten - A= FAINTING L, © 6, 6,Co, Y) 
where the function is given by the right side of equation 
(22). It is to be noted that the initial value of the 
relative phase difference, og, does not appear in the 


expression for the limiting value of (@-—-G). 


ete, 






3.  Daiscussion of Results. 

Equation (22) was solved numerically and the results 
are plotted in figure 1, which gives (#,_,—P) as a 
function of wavelength and drag coefficient. Using a 
surface geostrophic wind of 10 m/sec, the values of C,, in 
figure 1 correspond to frictional stresses of from zero to 
20 dynes/cm*. Cy=0.004 dyne sec/cm? corresponds to the 
value given by Sutton [4] for grass about 15 centimeters in 
height. The value of the static stability parameter, «= , 
of 2.90x107* cm? /dyne” corresponds to the value given by 
Gates [1] as the mean value over the United States in 
January. Those curves in figure 1 for which (®,- &) is 
zero when the drag coefficient is zero, are the curves for 
which term H in equation (11) is positive, i.e. these 
wavelengths are the stable waves since the requirement for 
instability in the frictionless case is that H be negative. 
The apparent amplification depicted in figure 2 for what 
should be stable waves with Cp)=0O can be explained by 
reference to equations (16) and (17) from which it can 
be seen that when C, and C_ are real numbers, the maximum 
possible value of B is given by B2=p* a BS st 2B_B The 
time variation of B for two stable waves, 3,000 km and 
12,000 km, with C,=0 is shown in figure 3. As is apparent 
from equation (17), B is a sinusoidal function of time, 
oscillating about the initial value. Figure 2 shows the 


effects of variations in the drag coefficient on the value 


10 


| ee 
a 7 = 
Bite 
-_ . 
_ _ — 
= > 
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bi 





of the amplitude for the wavelength band from 2,000 km to 
14,000 km. The values of B used to construct the figures 
are the values attained at the end of the first 36 hours 
for those wavelengths having (#,,—@® )>0, for C,=0, in 
figure 1. For the stable waves, due to the oscillatory 
nature of B mentioned above, the maximum value of B 
attained during the 36-hour period was used. For each wave- 
length the value of the initial phase difference between 
the temperature and pressure fields, oe, was such as to 
maximize the amplification in 36 hours (cc given to the 
nearest 45 degrees). As is obvious from figure 2, the 
effect of increasing the drag coefficient is to reduce the 
amplification for all wavelengths considered. This effect 
of friction is in agreement with the results presented by 
Holopainen [3] In comparing the results it should be 
noted that the ordinate in figure 2 is wavelength, whereas 
the ordinate in Holopainen's figure is a function of wave- 
length, static stability, and latitude. In figure 2 it 
can be seen that the wavelength for which the amplifi- 
cation is maximum decreases as the drag coefficient 
increases. 

Figures 4, 5, and 6 depict the effects of variations 
in thermal wind, static stability, and latitude on the 
maximum value of the amplitude for the same wavelengths as 
were used in figures 1 and 2, with Cy=0.004 dyne sec/cnt 


The central curves in figures 4, 5, and 6 are identical and 


Sh 





are the same as the middle curve in figure 2. From figure 4 
it is apparent that as the thermal wind increases the 
amplitude increases for all wavelengths, and that the wave- 
length for which the amplitude is a maximum also increases. 
Figure 5 shows that the effect of increasing the static 
stability is to decrease the amount of amplification for all 
wavelengths, except that for wavelengths greater than 
about 13,000 km the amplitude increases slightly. Also aso 
increases the wavelength of maximum amplification increases. 
The effects of moving the latitude northward are, as 
depicted in figure 6, to increase the amplification of the 
entire waveband, and to shift the wavelength of maximum 
amplitude toward lower values. It appears from the figures 
that menses the static stability and decreasing the 
latitude have essentially the same effects on the amplifi- 
cation, and that decreasing the thermal wind corresponds 
to increasing the drag coefficient in the resulting effect 
on the wave amplitude. One might conclude from the 
information presented in the figures that at low latitudes 
with a high value of the static stability parameter and a 
low thermal wind the pressure waves would not tend to 
develop, and that at high latitudes with low static 
stability and high thermal wind the tendency would be toward 
amplification of the pressure waves. 

Wiin-Nielsen [5) discusses the validity of applying the 


results of a linearized treatment of the dynamical equations 


12 








to obtain a limiting case for very long times. He justified 
the application on the basis that the initial trends were 
toward the situation given as the limiting case. With 
similar logic in mind, equations (17), (18), (20), and (21) 
were used to numerically investigate the values of the 
amplitude factors B and B,, and the relative phase difference 
(#%-¢d) for periods of time of less than 72 hours. It 
should be noted that while the limiting value of (@#-%) is 
independent of the initial phase difference, o¢, the short 
period variations of B, B+, and (@-¥) are very much 
dependent on the value of «. Figures 7 and 8 were 
constructed for a 5,000 km wavelength wave to depict the time 
variation of (¢#-4) and B, to show the dependence on e¢ of 
the time variation, and to show the fact that the terminal 


value of (@--¢Y) is independent of the value of ~. 


Because of the very similar behavior in the variations 
of B and B,- with time, only B is shown in the figures. 
Figures 7 and & may be compared with similar figures 
presented by Wiin-Nielsen [5], using the same wavelength, 
for the frictionless case. In comparing the two sets of 
curves consideration should be given to the differences 
in the values of the atmospheric parameters used in the 
calculations performed to construct the curves in figures 
7 and 8, and the values used by Wiin-Nielsen. One 
conclusion that can be drawn from the figures is that 
increasing the drag coefficient decreases the time required 
for (.- GP) to approach its limiting value. 
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